SLAC-TN-09-002 
BaBar Note # 327 
October 1996 



Noise in a Calorimeter 
Readout System 
Using Periodic Sampling 



Walter R. Innes 
SLAC National Accelerator Laboratory 



Abstract 

Fourier transform analysis of the calorimeter noise problem gives quantitative results on a) 
the time-height correlation, b) the effect of background on optimal shaping and on the ENC, 
c) sampling frequency requirements, and d) the relation between sampling frequency and the 
required quantization error. 
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Figure 1: Simplified photodiode-preamp circuit 

1 Introduction 

Noise in calorimeter readout electronics has been treated in some detail [Radeka, Haller, Dow]. 
However, since each author builds the details of their particular situation into their analysis, it is 
worthwhile to review the subject in light of our own experiment. Most previous studies assume 
that a single measurement will be made at a well determined time. Furthermore they assume that 
filtering preceding the digitization is limited to relatively simple analog devices. Since we have cho- 
sen to solve our lack of knowledge of the event time by recording periodic samples, applying digital 
signal processing techniques is natural. This also raises questions about the required sampling rate 
and the quantization error that don't occur in the usual scheme. 

Electronics noise and filters are usually analyzed using Fourier transforms [Ambrozny, Humphreys, 
Papoulis]. In the frequency domain, filtering, differentiation, integration, and time shifting are all 
represented by multiplications. Stochastic noise is most easily represented by its power spectral 
density (which is the Fourier transform of the more complicated auto-correlation function in the 
time domain). Many questions such as the signal to noise ratio for particular parameters and filter 
design can be handled entirely in the frequency domain. 

2 The Input Circuit Model 

Figure 1 shows a simplified version of the calorimeter input circuitry. Figure 2 shows a noise 
equivalent representation of the same circuit. The right hand portion represents the input FET. 
This is one of many possible equivalent representations [vanderZiel, page 29] [Ambrozny, page 137]. 
This particular one has a relatively simple connection between the physical processes which generate 
noise and the elements of the representation. The spectral densities of the noise generators will 
correspond to the those given in manufacturers specification sheets. 
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Photodiode JFET 



Figure 2: Model for the front end of the calorimeter readout. i s is the current due to the signal 
and is a function of time. Id is the spectral current noise generator corresponding to the shot 
noise of the photodiode. Ib is the equivalent current noise caused by background photons from lost 
particles. If is the FET input current noise, and V n is the FET input voltage noise. C s is the 
capacitance of the photodiode and Ci ss is the common source input capacitance of the FET. 



2.1 The photo-diode 

i s is a current generator corresponding to the signal generated by the interactions in the calorimeter. 
It will be represented by the function: 



.(*) = ~e-'/ T -lT(t). (1) 



t s is the decay time of the CsI(Tl) scintillation and H (t) is the Heaviside unit step function. t s is 
0.94 jus for CsI(Tl). The Fourier transform of i(t) is 

/OO 1 
i(t)e-3" t dt = — . (2) 
-oo 1 + JU}T S 

The spectral power density of the signal is: 

I* s (u)I s (u) = \ 2 . (3) 

For characteristic times represented as t x , I will use the notation that the associated radial velocity 
is uj x = 1/t x , the associated frequency is f x = uj x /(2tt), and the associated period is T x = l/f x . 
Using this notation we can write: 

J» = -=^- and I* s (u)I s (u;) = -J^. (4) 

Id is the current generator corresponding to the shot noise caused by the leakage current in the 
photodiode. This has a "white" (independent of to) spectrum with the value l\ = 2eli ea k (units of 
current squared per bandwidth). The order of magnitude for I\ ea k is nanoamperes. We'll assume 
our diode has a leakage of 4 nA. Since we average two diodes we can use the equivalent leakage of 
2nA for a single diode giving Id = 25fA/V Hz. 
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C s is the capacitance of the back-biased photodiode. A typical value is 80 pF, which we'll use 
in our estimates. 

A photodiode with these properties is the 1N2744. 

As much as possible, the discussion will be restricted to the electronic regime. However, some 
effects depend on the efficiency of the conversion of calorimeter shower energy into charge. When 
this number is required, I will use the value of 3000 photo-electrons per MeV in each diode. 

2.2 The FET 

Ci ss is the FET common source input capacitance. It can range from a few picofarads to tens 
of picofarads depending on the choice of FET. We will return to this choice later. In the circuit 
it functions in parallel with C s . Not shown in the figure since they are much less than C s are 
the feedback capacitance Cf and the calibration injection capacitance C C(l /. I will call the sum 
C s + C iss + Cf + C ca i simply C. 

If is a noise generator corresponding the input current noise of the FET. It is primarily due to 
the gate leakage current, although for very large gate areas the effect the thermal noise of the real 
part of the input impedance should be checked. The shot noise acts just like the diode shot noise 
and the two leakage currents may be added. In practice the FET gate leakage current is of the order 
of picoamperes and is negligible compared to the diode leakage current. Let l\ = 1% + ij « 1%. 

V n is a voltage noise generator representing the input voltage noise of the FET. This is pri- 
marily due the thermal noise in the FET channel. Many authors put a formula, which shows the 
dependence of this noise on the temperature and the forward transconductance, into their equa- 
tions [vanderZiel, page 75]. I prefer to keep V n as an explicit term. The value of V n is usually 
given on specification sheets and is also easy to measure. For junction FETs, V n is approximately 
independent of frequency over the range of interest to us. Typical values are 1 to 2 nV/ V Hz. 

In order to be more specific I'll take as an example the SANYO 2SK932 JFET. Ci ss is 20 pF 
and V n is « 0.7nV/\/Hz. Because we average the outputs two diode-FET-preamp combinations on 
each crystal, we divide by \f2 to get 0.5 nV/ V Hz. This is a good FET but not optimally matched 
to our problem. We'll return to this issue after we have developed the signal to noise formulas. 

2.3 Background 

lb is a current noise generator representing the sea of photons generated by particles lost from the 
beam. This noise resembles shot noise in that near impulsive elements arrive at random times. It 
differs in that the pulse shape is not a delta function but has the shape of the signal and there 
is a distribution of pulse areas instead of just the electron charge. 1 If we return to the derivation 
of the shot noise formula [Ambrozny, page 81], we see that the first difference can be treated by 
giving this noise the power spectrum of the signal (the Fourier transform of each pulse is the same 
as that of the signal except for phase and magnitude). The second is handled by noting that the 
shot noise formula is really a pulse area squared times a rate: t\ = 2e 2 fi ea k- e is the size of the 
pulse and fieak is the mean rate of their occurrence. Therefore we can get the low frequency limit 
of the noise spectral density of the background by doing the integral: 




(5) 



1 In radar theory, noise with these properties is called "clutter. 
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Figure 3: The amplitude spectra of the various noise sources expressed as input current densities. 
For the background, nominal conditions are assumed and no energy cutoff is used. Note that the 
signal has the same spectrum as the background noise. 

where B is the rate of background events per pulse area interval per time interval. The value of this 
integral can be inferred from the statistics of histograms such as those in BaBar Technical Design 
Report Figure 12-9 [BaBar]. Here I use more recent calculations of the background [Marsiske]. 
Assuming 3000 pe/ MeV/diode, the result for an average crystal is If, = 460 fA/VHz at nominal 
background. 

For convenience I'll use the symbol h to represent this low frequency limit and put the fre- 
quency dependence into the formulas explicitly. Since I n and V n are approximately independent of 
frequency, I can then treat /;>, I n , and V n as constants. 

Note that individual crystals can differ from the average by factors as large as 4. In particular, 
compared an average barrel crystal, an average endcap crystal has 5/7th the hit rate, twice the 
average energy deposit, and an 1^ twice as large. 

Figure 3 summarizes all the noise sources. 
2.3.1 Limitations 

The treatment of background as noise is a good description if the background signal after all 
processing contains significant contributions from many background events. If this is not the case, 
the data may fall into two classes, those showers little affected by background, and those with a 
large background contribution. The former will have a distribution of measured values consistent 
with electronic noise, while the latter may have a wider, highly skewed distribution. While the 
background as noise treatment may give a near optimal procedure for minimizing the rms error, 
it will widen the distribution for that class events unaffected by background. Depending on the 
physics objective, treating the background affected events as an inefficiency rather than events with 
errors may give the better result. 

How well does our case satisfy this many hit condition? In the background calculation in the 
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preceding paragraph, an average crystal had 0.11 hits/crystal/ fis. Assuming a window of 2 fis and 
a sum of 25 crystals, this implies 5.5 background hits contributing to the measurement during 
nominal conditions. This is enough to give reasonably Guassian behavior. However, much of lb 
is generated by the high energy tail of the background distribution. If we count down from high 
energies until get a mean number of hits of at least 2, the cut would be in the vicinity of 1 MeV. 
At 10 times nominal this cut would be close to 3 MeV. 

Another consideration is that an approach which tries to separate a background tail from other 
noise is hardly likely to succeed if the magnitude of the background contribution is of the same size 
as other errors. The size of the intrinsic error depends on the size of the signal. For a 100 MeV 
shower the expected error is ~ 2 MeV. Assuming that background contributions must be at least 
3 MeV before they can be separated from noise seems conservative. 

If only background hits with energies below 3 MeV are included, lb falls to 160 fA/V Hz. The 
background rate falls to 0.10 hits/crystal/ ^s. 



3 Theorems 



3.1 S/N 

The Radar Handbook [Skolnik, page 4-11] reviews the theory of finding a known signal in noise 
from the starting point of doing a least square fit to a finite set of measurements. Under fairly 
general conditions this is equivalent to the maximum likelihood fit and is the best that can be done. 
The treatment proceeds by taking the continuous infinite time limit which turns the usual matrix 
equations into integral equations. These are then solved by applying Fourier transforms. For a 
signal i(t, a±, a n ) with Fourier transform a±, a n ) the elements of the inverse of the error 
matrix for ai are given by: 

1_ r°° ff dljio) dTjuj) 
27r J~oo V V den da k 



1 r c 



/N(oj) duo, 



(6) 



where N is the double sided noise power spectral density. This matrix is known as the information 
matrix and also as the signal to noise ratio (S/N) matrix. 

In our case the signal has two parameters, amplitude and time offset: 



i(t) = ai (t - t ). 



This has the Fourier transform 



I{uj) = al (oj)e- jujt " . 
If we define Iq such that the actual value of a is 1 we find 

91 T< \ 



and 



(7) 
(8) 

(9) 



— = -jul(u). 



Thus the matrix of integrands in the S/N formula is 



1 



N(lo) 



(10) 



(11) 
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Because both |/(o;)| 2 and N(u>) are even in u>, the integral of the off-diagonal elements is 0. Since 
the matrix is diagonal, inverting it to get the error matrix V is trivial. If i(t) is normalized to unit 
area, the aa element of V is the square of the equivalent noise charge (ENC) and the to^o element is 
the square of product of the pulse height and the time error (ENTC). Since there is no correlation 
between these errors, a priori knowledge of the time does not improve the determination of the 
amplitude. The final formulas for the signal to noise are: 

'S_ 

and 

'S_ 

N 



1 Z 101 
2W-c 



N(u) 



du>, 



(12) 



1 

2^ 



N(uj) 



uj 2 doj 



(13) 



3.2 Matched filter 

This result for the error in the amplitude is the same that is given in many texts as the best that 
can be achieved using an optimal matched filter [Humphreys, page 65][Papoulis, page 135]. Assume 
there exists a filter such that the output peaks at t = 0. The square of that peak value is given by: 



o 2 (0) 



2 



(14) 



where F(u>) is the transfer function of the filter. The mean square noise after filtering is the same 
at all times and is given by 

< = ^j_ oo N ^)\ F ^)\ 2 du>. (15) 

The signal to noise is the ratio of o 2 (0) to n 2 . Multiplying the integrand of the numerator by 
y/N(u)/ ^jN{oj) and using Schwarz's inequality proves that the signal to noise is less than or equal 
to the S/N derived from the least square fit method (equation 12). Furthermore, inspection of the 
S/N equation before this manipulation shows that the equality can be achieved using a filter with 
a transfer function of 

F(co) oc Q^- . (16) 
TV (uj) 

From this I conclude that the best possible result can be achieved with a filter technique and that 
the required filter is readily calculated given knowledge of the shape of the signal and the spectrum 
of the noise. 

The impulse response to the optimal filter is 



1 r c 

'<'><**/-. 



and the shape of a noiseless signal after filtering (the convolution of the signal with impulse response) 
is 

(t)oc _L [°° JZP^du,. (18) 
v ' 2vr i-oo N(u) V ; 

A useful normalization for the filter can be achieved by dividing by the signal to noise ratio (equa- 
tion 12). This yields o(0) = 1. In the following discussion, F, f, and o are so normalized. 
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4 Application to our problem 



We have now collected all the necessary input information and the tools we need to design an 
optimal filter and calculate the signal to noise. Before tackling the full problem lets explore some 
special cases so that we can get an understanding of the effect of each input, can compare our 
results with previous work, and gain confidence in the method. Some of the formulas are followed 
by a bracketed reference to the source of the evaluation of the previous integral. Actually most 
of the integrals are straightforward (but sometimes tedious) to do using contour integration. This 
is a consequence of the fact that the signals, backgrounds, and analog filters all have exponential 
shapes in the time domain. More general signal shapes are less tractable. 



4.1 No background (I b = 0) and t s = 

This is the usual treatment when the "ballistic deficit" is ignored. i s (t) is a delta function, I s (u)) = 1, 
and 2N(lo) = I 2 + V 2 lu 2 C 2 . 

Before starting lets evaluate some useful constants related to the inputs so that we can use 
them to evaluate the results as we go: r n = V n C/I n = 2.0/is, co n = l/r n = 0.50Mradians/s, and 
f n = 80 kHz. 



S_ 

N 



ENC 2 



ENTC 2 
F{w) 



Ttl 2 J-oo OJ 2 + OJ 2 



doj 



- s~ 








tt 





I 2 



In T n — InVnC 



oo tol + ^ 2 



duj 



oo 



(19) 

(20) 

(21) 

(22) 

(23) 
(24) 

(25) 



/(*) 
o(t) 



- p -\t\/Tn 



- p -\t\/T„ 



(26) 

(27) 



Substituting our sample values, we find that ENC = 223 e, and the corner frequency of the 
optimal filter is 80 kHz. 
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4.2 No background (I b = 0) and r s > 

This is the case usually treated, ballistic deficit included. 



hi") 

2N(u) 



S_ 

N 



ENC 2 

'S_' 

N\tt 

ENTC 2 
F{u) 

/(*) 



o(t) 



U - JUJ S 

r2 I T7"2, ,2/~y2 



(i 2 + v 2 ^c 2 ) = i 2 (^ + ^ n )M 
i 



< r 
2 

n J—< 



VT/2 y-ooC^ +^(0,2+^2) 



do; 



[Dwight 856.31] 



co 2 u 2 >' x 



'1 f 
2 / 

T). J—' 



0J 



7T/2 y_ 00 ( W 2 +w 2 )(ti; 2 +ti; 2 : 

^r n r s (r n + r s ) = ENC 2 r n r s 

2juj n (uj n + u s ) 
(uj+juj s )(u} 2 +u>l) 



da; 



[Maple] 



(28) 
(29) 
(30) 

(31) 
(32) 
(33) 
(34) 
(35) 

(36) 



-oo [U 2 + U 2 S ){UJ 2 + U) 2 

U ^ n + ^ ) (e-W* - ^e-l'l/- - siguum(t) U~W* - e*)) (37) 



( e -\t\/rn +H{ _ t) J^( e -\t\/r n 
V T n - T s V 



-\t\/Ts 



)) 



7T 



(0,2+^2)^2+^2) 



djjj 



(r n - t s ) 



-\t\/rs 



) 



(38) 
(39) 
(40) 



Substituting our sample values, we find that ENC = 239 e, ENTC = 327/Use (1.1ns for a 
100 MeV deposit), and the corner frequency of the optimal filter is still 80kHz. The ENC corre- 
sponds to and equivalent noise energy (ENE) of 80 keV. 
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4.3 I h > 0, V n = 0, and I n = 



This special case of no electronic noise is explored to give us confidence that the filter technique 
does separate the desired signal from the background noise. 



2N(lo) 
\I I 2 



N 



U - ju> s 

, , 2 T 2 
UJ s I b 

UJ 2 + LO 2 



S_ 

N 



= oo , ENC 2 = 

= oo , ENCT 2 = 
1 f°° 

f(t) OC - / (Us-jLUy^dw 



S_ 

N 



oc 2(5(t)-T s 5'(t)) 



o(t) oc 



7T 



e jujt duj oc 2<5(t) 



(41) 

(42) 

(43) 

(44) 

(45) 

(46) 
(47) 
(48) 



Because the signal and (background) noise have the same spectrum, the signal to noise ratio is 
independent of frequency. This results in infinite signal to noise integrals and no error on the pulse 
height and time determinations. In the absence of other noise, the matched filter turns both the 
signal and background events into delta functions which are always separable. 

4.4 The general case 

Now we will treat our general case, with all noise and width processes active. The signal is as for 
the previous two cases. The noise is the sum of the noise in those cases. 



|/ S HI 2 

2N(u) 



U - juj s 



UJ 2 + LO 2 



T 2 i , 2 T 2 

n (, , 2 4. , ,2^ i ^sJ-b 

-tK+ w ) + 



(49) 
(50) 
(51) 
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s_ 

N 



S_ 
N 



UJ 2 U) 2 < x 

s ■ n 



fit) 



(t) 



1 



b 2 = 



VT/2 W 4 + (a; 2 + u 2 )(jJ 2 + W 2 W 2(1 + 72^2) 

1 + if /I 2 and w 6 = w n 6 



duj 



(Il + lD^Us + ^l+uj) 



[Dwight 857.11] 



£iVC 2 = {I 2 n + lt)^Wsp + T b /T s + T s /(b 2 T b ) 



- s~ 








tt 


irl 2 J_ 



7T/2 W 4 + (W 2 + w2> 2 + W 2 W 2 (1 + /2/J2) 



duj 



(52) 

(53) 
(54) 

(55) 
(56) 



ENTC 2 = I 2 n T n T s ^T 2 + 7-2 + 2bT n T s ) = ENC 2 T b T s [Prudnikov 2.2.10 #4 page 313] (57) 
F{u) = ENC 2 

fit) 



2uj s uj 2 1 (uj s +juj) 

I 2 (^ + (u 2 +U 2 n )u 2 +U 2 U 2 n b 2 ) 

2 "•" {u s + juj)e jujt 



ENC 2 **** ' ' X 



uj s cos{ujt) — to sin(u;t) 



duj 



-duo 



u+ = 



(58) 

(59) 

(60) 
(61) 



ENC 2 ^ / 
ttI 2 J- 

TTlf {l ( W 4 + ( W 2 + w2> 2 + ^2^) 

v / (^ s + K 2 + ^ 2 t )/2)/2 
^{u> b u s -(u 2 + u 2 )/2)/2 (62) 
ENC 2 ^ e -+l'l f J- cos(u;_i) + -L sin(w_|t|) - sin(c^t)") (63) 

[Prudnikov 2.5.10 #15 and#17 page 397] 



-\t\/r+ 



ENC 2 



-\t\/r+ 



( cos(i/r_) + ( — signum(t) — ) sin(i/r. 

V \T+ 7"b/ 



,2, ,2 



7T/2 



(W 4 + K 2 + ^2)^2 + ^2^2) 

cos(i/r_) + ^sin(|t|/T_)) 



duj 



(64) 
(65) 
(66) 



The filter output for a noiseless signal pulse is shown in Figure 4. 
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Figure 4: A noiseless signal after passing through the optimal filter. The peak gives the estimate 
for the signal charge and time. In this case t=0 and the charge is unity. 

4.5 I n = 0, u s « ujh 

The results in the previous section are a bit complex, but can be simplified for our case if we 
note that If> is more than 10 times even at nominal background levels. In addition, uj[> is 
more than five times uj s . Setting I n = 0, using the approximation uj^ » uj s , and defining 
uj s b = rs ijj^ k, lo + , we find: 



-ju s 



(67) 



LO - JLO s 



h(u)\ 2 = 



UJ 2 + LOg 



(68) 



2N(lo) 



.) 



(69) 



I s {u)\ 2 /N{u) 




(70) 




(71) 




(72) 
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Figure 5: The optimal filter in the frequency domain. 
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Figure 6: The optimal filter in the time domain. 



F{u) 



fit) 



o(t) 



ju b UJ sb (uJ - juJ s 



\t\/r BI 



cos(t/T sb ) +sin(|t|/r s6 ) - \2T s /r b sin(t/r. 



r *6. 



l*l/ T " (cos(t/r sb ) + S m(\t\/r sb )) 



(73) 

(74) 
(75) 



The optimal filter is shown in the frequency and time domains in Figures 5 and 6. 

Now for some values: r b = V n C/I b = 0.11 /is, the filter corner yJr s T b = 0.32 /j,s, ENC = 1930 e, 
and ENTC = 620 /us e or 2.1 ns for a 100 MeV deposit. The frequency corner of the optimum filter 
is 0.50 MHz. 

For our worst case of times 10 x nominal background, goes up x3.1, r b goes down x3.1, and 
y'ryrr, goes down xl.78 to 0.17 fis. ENC = 4460 e, ENTC = 760 fise, and the corner frequency 
goes to 0.90 MHz. If the 3 MeV cut were used to calculate I b , I b would be lower by a factor of 2.4 
and the values for x 10 would be about the same as the uncut nominal values. 
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5 Picking the best FET 



We see that the ENC depends on the FET properties only in the combination V n (Ci ss +C s ). The 
smaller this value the better. For a given FET design, these parameters can be varied by changing 
the gate area A (or equivalently connecting FETs in parallel) . V n decreases as 1/ \[A~ while Ci ss is 
proportional to A. V n C is a minimum with respect to A when Ci ss = C s . This implies that our 
sample FET should be scaled to 4 times its original gate area. Making this change would decrease 
V n C by a factor of 1.25. Since ENC goes as the fourth root of V n C, ENC would decrease only 
5%. On the other hand, ENTC would decrease by 16%. 

6 Sampling Requirements 

So far we have assumed that the signals are ideally and continuously measured. In fact they will 
be sampled with limited precision. What restrictions on the sampling rate, length, and accuracy 
will prevent the loss of information? 

6.1 Sampling frequency 

The sampling frequency must be high enough to capture all frequencies with useful S/N. It must 
also be high enough to avoid misrepresenting noise at higher frequencies as noise in the signal region 
(aliasing). On the other hand the sampling frequency should be no higher than necessary. 

What is the highest useful frequency? We note that S/N is approximately constant up to 
to = uj s b and then decreases as 1/uj 4 . If we throw away all information above w = uj c the fraction of 
information lost is approximately (3/4)(u; s 5/u; c ) 3 . For a 20% loss of information this implies that 
to c > 1.5lu s }), and that f c > 0.6 MHz. The fractional loss in time signal noise is somewhat larger, 
being given approximately by (3/4)(u; s b/u; c ). Thus our choice based on the amplitude error will lose 
us half of our time information and increase our time error by 25%. This is probably acceptable. 
The exact formula for the fractional loss of amplitude information is 

1 _ I 1 1 i og I ^j+2^c + 2o4 
■k\2 2io sb LO c + 2u? sh 

Assume that we have one stage of near ideal integration before the digitizer (the charge sensitive 
preamp) followed by a near ideal differentiator. This pair gives a band pass filter with a pass band 
from a low frequency to f c . The fall in the stop band is 1/uj. Further assume that there is an 
additional two pole low pass filter with a corner at f c (drop off like (u c /ll>) 2 for uj > uj c ) between 
the preamp and the digitizer. If we sample at frequency f d , noise at frequencies above fd/2 will 
appear at f b s = fd~ f ■ We are interested in studying the contribution to N from V 2 which drops 
off as (ujc/w) 4 after all the analog filters. We have sufficient protection against aliasing if V 2 at 
/d/2 - f c « Vn^Jc- This is satisfied if (f c /(f d /2 - / c )) 4 < 4. This implies f d > 3.4/ c , which 
is 3.4(1. 5)/(2TTy/T s n) = This is 2.5 MHz at nominal background and 4.7MHz at xlO 

nominal. 3.7 MHz is adequate adequate for nominal backgrounds and marginal for xlO. 

6.2 Sample length 

The number of samples required is related to the lowest frequency of interest. So far we have 
assumed that all integrals go from zero frequency. Pile up makes this impractical as does the fact 



+ tan 



2L0 c L0 s b 



(2c 




(76) 
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that we can't use all past history and wait forever to get an answer. 

Lets look at the pile up problem first. Allowing for variations from crystal to crystal, and warm 
diodes, we should plan for at least 20 nA of total current in a diode. If the offset from this current 
is to stay within 10% of full scale, the product of this current and the integration time must be 
less than 10% of the charge of a 14GeV shower, i.e., < 0.5 pC. This limits the integration time to 
200 us. This sets the minimum value for the frequency corner of the first integrator. 

This gives 740 samples as the maximum useful number at 3.7 MHz. This is a large number. 
What is the smallest number of samples we can use without significant loss of information? For 
the amplitude parameter, the S/N is flat up to the corner frequency and then drops rapidly. For 
our lowest noise case (ij, = during source calibration for example) , we have a corner frequency of 
80 kHz. If we permit a 20% loss of information, we can cut off the integral at the low frequency 
side at 13 kHz. This is less than three times the previous integration limit leaving us little choice 
and many samples. The problem here is that the sample rate is much higher than needed for this 
case. Once filtering has been performed, the sample set could be decimated, but most of the work 
is done by then. In the nominal background case the corner frequency is 400 kHz, and the integral 
could be cut off at 80 kHz for a sample length of 12.5 (is or 46 samples. This would be the normal 
operating mode. If the sampling frequency were optimized for this nominal background level, the 
number of samples would be less than 25. 

6.3 Sampling summary 

We can summarize the conclusions of the last two sections by noting that the number samples that 
must be dealt with is proportional to the ratio of the highest and lowest frequencies of used. The 
proportionality constant lies somewhere between 2.5 and 5 depending of the details of the system 
and how "highest" is defined. If the signal to noise ratio is approximately constant from low 
frequencies up to the highest frequency of interest then the fraction of the available information 
retained (the efficiency) of the system is given by 1 — f lowest I f highest- An interesting but not 
necessarily relevant observation: the maximum information per sample is obtained when fi owes t is 
one half of f highest and the efficiency is 50%. 

7 Quantization Error 

7.1 With no net integration or differentiation 

With the same assumptions about filtering as in the previous section we can estimate the require- 
ments on the quantization error. We assume there is a low pass filter such that there is no aliasing, 
no loss of information, and there are offsetting integrators and differentiators such that there is no 
net effect in the pass band. Under such conditions, quantization error may be referred back to the 
input where it appears as another current noise source. 

According to the sampling theorem [Papoulis, page 141], if the anti-aliasing conditions are met, 
the original signal may be reconstructed from the samples by the interpolation formula 



sm((u d /2)(t -nT d )) 
(co d /2)(t - nT d ) 
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where is the sampling period, and uj^ is related to it in the usual way. An error in a sample 
can be represented as a function proportional to one element of this sum. The Fourier transform 
of such a function is flat out to fa/2. The amplitude of this noise pulse is drawn from a square 
distribution whose width is given by the bit resolution (or the effective bit resolution for a non-ideal 
ADC). The rms area of the pulse is then T^i^ / \/l2. The rate of pulses with effectively random 
phase is fd = l/Tj, and the noise spectral density is flat with a value of I 2 = i^ t /(12/^). This is to 
be added to the FET input current noise. We should compare this to I 2 for very quiet conditions, 
and to i| for nominal conditions. 

First lets treat the quiet case for which ENC 2 goes as \fT^. If we wish to restrict our error 
increase to less than 5% due to this source we require that I 2 < I 2 /5. Therefore i^t should be less 
than y/2fd I n - For the example diode and our sampling rate, this is %n < 68 pA. The current at 
the signal peak is q/r s if there is no analog shaping. For a triple RC filter with a time constant 
of 0.25 ^s, this is reduced by a factor of 2. Assuming that 12GeV somehow gets into one crystal, 
the maximum peak current would be 3 fiA. This implies a dynamic range requirement of 43,000 
or 15.4 bits. Since the signal is not calibrated, we need to add a bit for gain variations suggesting 
that 16.5 effective bits of dynamic range are required during quiet conditions. Note that better 
light collection imposes a greater dynamic range requirement. 

During nominal conditions, ENC 2 goes as The 5% error increase condition implies that 

I 2 < I 2 /7, and that iut be equivalent to less than y/TTTJah- For nominal background conditions 
and our sampling rate this is i^t < 710 pA. The dynamic range requirement is 4,200 or 12 bits. 
With an extra bit for calibration differences, this is 13 bits. 

In practice, many crystals are summed to measure a shower. Both the electronics and back- 
ground noise increase as y/m, where m is the number of crystals included. The dynamic range 
required decreases as \fm. For a 9 crystal sum and nominal background conditions, the number of 
bits required is 10.5. 

7.2 With one net integration 

If we assume that instead of offsetting integrators and differentiators, we have only one integrator 
and then the low pass filter, the quantization error behaves as a voltage noise when referred to the 
input. The rms area of the error pulse is T^vutl 'V^2, and the noise spectral density is flat with 
a value of V 2 = v 2 it /(12fd). This is to be added to the FET input voltage noise V 2 . Recall that 
in our case ENC 2 goes as y/V 2 in the quiet case and as \/V 2 in the nominal case. If we wish 
to restrict our error increase to less than 5% due to this source we require that V 2 < V 2 /10 in 
the quiet case and V 2 < V 2 /5 in the nominal case. This implies that vut be equivalent to less 
than y/1.2fd V n and \/2Afd V n , respectively. For the example FET and our sampling rate this is 
Vfnt < 1.05 /iV and vut < 1.5 /iV, which given the source capacitance corresponds to a charges of 
< 660 and < 940 e, and to a shower energies of < 0.22 and < 0.31 MeV. Assuming that 12 GeV 
somehow gets into one crystal, the dynamic ranges required are 55,000 and 42,000, or 15.7 bits for 
quiet conditions and 15.4 bits for nominal conditions. Since the signal is not calibrated, we need to 
add a bit for gain variations suggesting that 17 effective bits of dynamic range may be sufficient. 

All of the above dynamic range calculations address the sufficient conditions for which the 
quantization error will not contribute to the electronic (and background) noise. They do not treat 
whether or not this dynamic range is necessary given the inherent energy fluctuations of the shower 
process. 
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The considerations for multi-range digitizing have been addressed by Al Eisner and Gunther 
Haller and remain unchanged. For larger pulses, sources of error other than electronics dominate. 



8 Filter Implementation and Interpolation 

There is still the problem of signal extraction. Because it is non-causal, a matched filter cannot 
be implemented as a simple analog filter. A procedure for directly calculating the coefficients of a 
tapped delay line approximation to the optimum matched filter is described by Papoulis [Papoulis, 
page 327]. The result is exactly parallel to the usual least square fit formulas. The autocorrelation 
function, rif, is the inverse Fourier transform of the noise power spectral density after the analog 
filter. The elements of r, the error matrix for the samples, are given by nj(\l — k\T). The expected 
values a for the samples is given by i s f((l — m)T + 1), where i s f is the expected signal after analog 
filtering and t is the time offset of the signal from sample m. The coefficients of the tapped delay 
line filter are given by r _1 a. The optimally filtered estimate for the signal at time t is given by 
o(t) = dr *a, where d is the set of samples. For the case of the triple RC analog filter, the 
coefficients each have four terms with different t dependencies. The estimate for o(t) also has four 
such terms and its derivative with respect to t has three. Calculating o(t) takes four multiply and 
adds for each sample used. Finding the root of the expression for the derivative with Newton's 
method should not take many iterations. Each iteration involves calculating one exponential and 
approximately a dozen multiply and adds. The presence of multiple peaks in the time window or 
the absence of any peaks would complicate this last step. A possibility is to not try to find the 
peak, but to report the coefficients of the four terms instead. 

Knowing the event time would reduce the time to get an estimate for o(t) by a factor of four, 
and there would no need to search for the maximum. 

9 Less Than the Best 

The previous sections deal largely with optimum solutions, albeit with some practical limitations. 
This section will examine the loss of information if other that optimal filtering is used. It may be 
less than optimal in that it is not matched, e.g., an all analog filter, or in that the pass band is 
optimized for another condition. 

Let us look at the case of an all analog filter. In this case all filtering is done prior to digitizing 
and the only digital processing performed is interpolation. I will take the BaBar TDR design 
(described in subsection 6.1). This has a charge sensitive preamp followed by a CR and then two 
RC filters. Together these give the equivalent of three RC filters, all with the corner frequency f c . 
This has the transfer function 



The time dependence of the output signal is given by the inverse Fourier transform of their product. 
This integral can be done by contour integration. If the contour is completed by a semi-circle at 
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Figure 7: Raw input current signal and the 
signal after triple RC niters with time con- 
stants of 0.25 [is and 2.0 [is. 



Figure 8: The peak signal current after a 
triple RC filter vs. the filter time constant 



infinity around the upper half plane (for t > 0), the integral can be deduced from the residues at 
the poles jio c and ju s . The result is 

*> ^ lA { " <«■ " ^ + ') " e """} ■ (80> 

This is shown for several time constants in figure 7. The peak value vs. time constant is shown in 
figure 8. Since there are no poles in the lower half plane, the integral is for t < 0. 

The mean square noise is given by the inverse Fourier transform of the noise evaluated at t = 0, 
i.e., the autocorrelation for no time difference. The noise before filtering is given by equation 51. 
The mean square noise is 



1 



oo 
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n=- N{uj)\F{uj)\ 2 duo. (81) 
vr Jo 

The contribution of the electronic noise is 

Tlelec = Uc(VnC 2 iu 2 c + 3I 2 )/32. (82) 

The peak of \/n e i ec /o 2 (t) gives the ENC for the electronic noise after analog filtering. Appropriate 
scaling gives the ENE shown in Figure 9. The background noise contribution is 



(8uj 2 + 9uj s u; c + 3lo 2 )uj c uj s I 2 
32(u c + uj s y 



nback = j l°3J 



n is given by their sum. The signal to noise ratio for a single sample taken at time t is given 
by o 2 (t)/n. If we divide this by the optimum signal to noise given by equation 55 we get the 
"efficiency" for this single sample versus the sample time. Figure 10 shows this function for the 
BaBar calorimeter during nominal operating conditions for several values of the corner frequency. 
The highest efficiency achieved is 68% for a shaping time constant of 0.2 fis. For low background 
conditions (Figure 11) the peak efficiency is 80% at a shaping time of 1.5 [is. At ten times nominal 
background, (Figure 12) the efficiency drops to 50% at a shaping time of 0.12 [is. The loss of 
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Figure 9: The Equivalent Noise Energy (elec- 
tronic only) of a single sample at the peak af- 
ter a triple RC filter versus the radial velocity 
of the filter corner (1/RC). 



Figure 10: The efficiency vs. 1/RC for nomi- 
nal background conditions. 



efficiency for nominal and better background conditions is not significant. The loss at during high 
background conditions hits just when noise is the causing the most problems. Matched filters utilize 
knowledge of the signal shape. In low background conditions, the optimum shaping time is longer 
than the signal. The shape of the signal is not seen and matched filters offer little advantage. The 
situation and conclusions are reversed in high background conditions. 

The cost of using an analog filter optimized for conditions other than those encountered is more 
dramatic. If the corner frequency is set to the 0.2 fj,s which is optimal for nominal conditions, but 
the backgrounds are actually at 10 x nominal levels, the efficiency is 44%. If the corner frequency 
were optimal for low background conditions, the efficiency at 10 x nominal background would be 
12%. 

10 Summary of Significant Findings 

Since the important conclusions may have gotten lost in the derivations, I'll review them here. 

• Machine background can be treated as a current noise with the spectrum of the signal. 

• The best possible amplitude and time measurements can be obtained with a matched filter. 

• The errors in the amplitude and time are not correlated, therefore knowing the time does 
not improve the ultimate precision of the amplitude determination, although such knowledge 
may reduce the computation required. 

• The signal to noise methodology used here reproduces previous results. 

• The background current noise dominates the photodiode shot noise even at nominal back- 
grounds and in the quietest part of the calorimeter. 



18 




4 6 
1/RC (Mrad/s) 



10 




1/RC (Mrad/s) 



Figure 11: The efficiency vs. 1/RC for low 
background conditions. 



Figure 12: The efficiency vs. 1/RC for 10 x 
nominal background conditions. 



For nominal background where Ib^> I n - 



ENC 2 « ltJ2T a n + r* = j2T s llV n C + I 2 b V^ 



ENTC 2 « ENC 2 T s r h = J2I b rfV^C 3 + r|F n 4 C 4 



• The sampling frequency of 3.7 MHz is adequate at nominal backgraounds, but it is marginal 
for lOx nominal background levels. 

• At this frequency, 64 samples will have to be processed under nominal conditions. 

• The step of an ADC bit (on the most sensitive range) should be set so that there are at least 
13 bits to full scale (12GeV) during nominal background. For quiet conditions and with a 
differentiating stage in the analog filter, 16.5 bits may be required. 

• Simple analog filters can achieve near optimal results in low to moderate background condi- 
tions if the shaping time is optimized for the actual conditions. 
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